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We present the leptonic decay constants fj)^ and /£,+ computed on the MILC collaboration's 2+1 
flavor asqtad gauge ensembles. We use clover heavy quarks with the Fermilab interpretation and 
improved staggered light quarks. The simultaneous chiral and continuum extrapolation, which 
^—1 determines both decay constants, includes partially-quenched lattice results at lattice spacings 

On a w 0.09, 0.12 and 0.15 fm. We have made several recent improvements in our analysis: a) we 

^P, include terms in the fit describing leading order heavy-quark discretization effects, b) we have 

^ adopted a more precise input ri value consistent with our other D and B meson studies, c) we 

k> have retuned the input bare charm masses based upon the new ri . Our preUminary results are 

/o^ =260±10MeVand/o+ =217±10MeV. 
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1. Introduction 

We report on progress in the Fermilab Lattice and MILC Collaboration calculation of the 
D meson decay constants. This work is a continuation of the program that predicted the decay 
constants: /0+ = 201(3)(17) and fo, = 249(3)(16)MeV [1], in good agreement with the CLEO-c 
value of fj)+ = 205.8 ±8.5 ±2.5 [2, 3]. We have since extended this calculation to two additional 
ensembles at our finest lattice spacing a 0.09 fm and we have replaced a limited set of very coarse 
(a « 0.18 fm) ensembles with higher statistic ensembles at a somewhat finer spacing a ^0.15 fm. 
In our last update [4] we reported: /£>+ = 207(11) and fo^ = 249(11), where fo^ is about 0.6c7 
lower than the recent experimental average [5]. The value of fo, remains an pressing issue given 
that experimental average is about 2. la higher than the most precise lattice result from the HPQCD 
collaboration [6] . The apparent tension between experiment and lattice predictions has motivated 
suggestions of physics beyond the Standard Model [7]. 

Smaller statistical uncertainties and better control of systematic effects are key to resolving 
the puzzle. In this report, we have doubled statistics on the most chiral of the a 0.09 fm 
lattices; otherwise, statistics have not changed. A new generation of calculations, now underway, 
aims to increase statistics by a factor of four overall. Our progress includes: a) a better method 
of accounting for heavy-quark discretization effects b) a more precise input value for the scale 
parameter ri, consistent with our other heavy quark studies and c) more precisely tuned input 
charm kappa values. 

2. Staggered chiral perturbation theory for heavy-light mesons 

We use the asqtad improved staggered action for both sea and light valence quarks. Leading 
discretization effects split the light pseudoscalar meson masses, 

^ab,^ = {f^a + mb)n+ a^A^ , (2. 1) 

where there are sixteen tastes in representations ^ = P,A,T,V,I. 

Staggered chiral perturbation theory for heavy-light mesons accounts for such taste breaking 
effects [8]. At NLO in the chiral expansion, for 2 ± 1 flavors, and at leading order in the heavy 
quark expansion, 

= <J>o [1 + AfH{mg,mi,mh) + pH{mg,mi,mh) +CLaya^] , (2.2) 

where (pH.^ = fuq -v/"^ ^^^^ fiiq is the decay constant of a heavy meson Hq consisting of a heavy 
quark and a light quark of mass mq. The heavier sea quark has mass m/, and the two degenerate light 
sea quarks have mass m/. The <pH,,, in general, are partially quenched: niq ^ mi and niq ^ mf,. The 
chiral logarithm terms, Afy, are a dependent as a consequence of the mass splittings in Eq. (2.1) 
as well as from "hairpin" terms proportional to the low energy constants a^5^ and a^d^,. The a 
dependence of the analytic terms, pn, ensures that is unchanged by a change in the chiral scale, 
A;^, of the logarithms. The expression in Eq. (2.2) is used in our combined chiral and continuum 
extrapolations. In practice, we add the NNLO analytic terms to the fit function in order to extend 
the fit up to niq nis and extract fo^. Priors for the parameters a^5^ and a^d^, as well as values of 
the physical light quark masses are obtained from the MILC analysis of fji and fa [9]. 
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3. Discretization effects from clover heavy quarks 

We use tadpole-improved clover charm quarks. At leading order, discretization errors are a 
combination of ^(a^A^g) and &{aaKuQ) effects where a is the QCD coupling and A//g is the 
scale in the heavy quark expansion. Our past studies have estimated heavy quark discretization 
effects using such power counting arguments to bound the error at the smallest lattice spacing, 
taking Auq ~ 700 MeV. This rather crude method does not effectively use the data to guide the 
error estimate. 

This study introduces a new procedure: the leading order heavy quark discretization errors are 
modeled to leading order as part of the combined chiral and continuum extrapolation. At tree-level, 
discretization effects arise from both the quark action and the (improved) current [10, 11, 12]. We 
add five extra terms to Equation 2.2: 

[a^AHQ{cEfE{amQ)+cxfx{amQ) + cyfYiamQ)} + avaAHQ{cBfB{amQ) + C3f3{amQ)}] 

(3.1) 

The coefficients ce, cx, cy, cb and C3 are additional parameters determined in the fit while the /, 
ai^e (smooth) functions of the heavy quark mass, aniQ, known at tree level. We introduce priors for 
the coefficients constraining them to be ^(1) while setting Khq = 700 MeV and nic ~ 1.2 GeV. 
Currently the data are too noisy and the shapes of the functions /, are too similar for the fit to prefer 
a particular A//g. Including the heavy-quark discretization terms increases the decay constants 
by a few MeV and increases the eiTor from ~ 1.8% to ~ 3.8%. The larger enm now includes 
the residual heavy-quark discretization uncertainty in addition to residual light-quark discretization 
effects (encoded in Eq. (2.2)) as well as statistical errors. We quote a combined uncertainty from 
all the three sources of error. 

4. Lattice spacing determination from r\ 

The distance r\ is a property of the QCD potential between heavy quarks. The ratio r\/a, for 
lattice spacing a, has been computed for all of the MILC gauge ensembles. At intermediate stages 
of the decay constant analysis quantities are converted from lattice units to r\ units using r\/a. 
The value of r\ must then be input in order to convert results to physical units. The r\ value is 
also an input to the process of determining other quantities such as the bare charm quark masses as 
discussed in the next section. 

Figure 1 depicts several r\ determinations. The first two determinations historically (circa 
2004-2005) are labeled "HPQCD T{2S-IS)" [13] and "MILC T(2S-1S)" [14]. They are both 
based on the same analysis of the T spectrum by the HPQCD Collaboration using a subset of the 
current MILC ensembles. The two determinations differ mainly in the details of the continuum 
extrapolation. The MILC Collaboration is also able to infer a value of ri based on the value of fji 
they find in their analysis of the light mesons. Recent light-meson analyses include results from 
finer lattice spacings than the earlier T spectrum study and the resulting ri values are known to 
better precision. The figure shows the result of two recent analyses labeled 'MILC 2007" [15] 
and "MILC 2009" [9]. The (preliminary) 2009 result agrees at the 0.9a level with the MILC 
T value but differs from the HPQCD T value at the 1.8 a level. As these proceedings were being 
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r1 [fm] 

MILC Y(2S-1S)- I I I 

HPQCD Y(2S-1 S) - I I 1 

MILC fpi 2007- I I I 

MILC fpi 2009- I I I 

H PQCD 2009 - I I I 

I I I I I 

0.305 0.310 0.315 0.320 0.325 

r1 

Figure 1: Values of scale parameter ri in fermi units. The "HPQCD Y(25-15)" value uses the HPQCD 
collaboration T spectrum results to set the physical value [13]. The "MILC T(25'-15')" value derives from 
essentially the same spectrum analysis [14]. MILC determines ri more precisely from their calculation of 
f^: "MILC /;r 2007" [15] and "MILC 2009" [9]. In a very recent update, "HPQCD 2009", several 
physical quantities, including recent T results, are used as inputs [16]. 



ri [fm]: 


0.3108 


0.318 


a K run 


K tune 5<ps % <ps 


K tune Sips % <ps 


0.09 0.127 
0.12 0.122 
0.15 0.122 


0.1272 -0.0043 -0.56 
0.1222 -0.0036 -0.50 
0.1222 -0.0031 -0.42 


0.1267 +0.0065 +0.84 
0.1215 +0.0091 +1.26 
0.1213 +0.0108 +1.47 


[MeV] 


-1.8 


+ 1.3 



Table 1: Tuning of k charm at the three lattice spacings for two choices of ri. The shift 5(ps is the change 
in (j) at the strange quark mass when K changes from the run value to tuned K value. The corresponding 
change in extrapolated /^^ is 5fo^. In each case, all other extrapolation inputs are fixed to their appropriate 
(ri dependent) values. 

prepared, HPQCD published a new value for r\ [16], labeled "HPQCD 2009" in the figure, in much 
better agreement with IVIILC's recent ri values. 

In this study, we use the IVIILC ri determinations from /^r to set the physical scale. Our 
central value for ri (the 2007 value) was also used in our studies of the semileptonic decays on the 
same lattices [17, 18]. The range of the 2009 IVIILC ri determination is used to set a symmetric 
uncertainty around the central value. Hence, we take ri = 0.3108 ±0.0022. Our previous decay 
constant studies used the IVIILC T value: ri = 0.318 ±0.007 as an input which is about one a 
higher. 

5. Retuning kappa charm 

We determine the value of K for the charm quark by requiring that the spin-averaged kinetic 
masses of the lattice pseudoscalar and vector mesons made from a heavy clover quark and strange 
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0Dj 


0Z)+ 








3 4 


1 3 


chiral extrapolation 


2.0 


2.5 


0.8 


inputs ri, m^, and m„ 


0.7 


0.7 


0.3 


input 


1.0 


1.2 


0.2 


Z^' and 


1.4 


1.4 





higher-order 


0.3 


0.3 


0.2 


finite volume 


0.2 


0.6 


0.6 


total 


3.8 


4.7 


1.7 



Table 2: Uncertainties as a percentage of and the ratio. The total combines all of the errors in quadrature. 

asqtad valence quark equal the spin-averaged Dg meson mass. The tuning depends upon r\ in the 
conversion between lattice and physical masses. 

In the past year, we have conducted new kappa-tuning runs with at least four times the statistics 
of our older tunings. At each lattice spacing, we simulated mesons for three values of K around 
charm and three light-quark masses around strange allowing us to retune K for a given r\ . 

Table 1 shows preliminary tunings for K charm based upon the two input values: ri = 0.3 108 fm 
(present value) and r\ = 0.318 fm (past studies). For each ri, the (preliminary) tuned kappa and the 
corresponding change 5<^s = (j)s{^ tune) — run) is listed by lattice spacing. The "run" kappa 
values are those used in the decay constant simulations. We adjust each (j)q point by 50., prior to 
the chiral extrapolation to correct for the mistuning of kappa. The bottom row of the table shows 
the resulting change in fo^ . The opposite signs of the differences show that keeping kappa tuned 
partly compensates the change in r\. We find that changing ri from 0.318 fm to 0.3108 fm while 
keeping kappa charm tuned increases fo, by about 4.2 MeV. 

6. The chiral and continuum extrapolation, results and uncertainty budget 

We fit (pDg results from lattice simulations on eleven asqtad MILC ensembles [14] at the three 
lattice spacings: a w 0.09, 0.12 and 0.15 fm. Our valence quark masses are in the range O.lm,, < 
niq < nis. Since our last report, we have doubled the statistics at the most chiral of the a ^ 0.09 
ensembles. The 3x4 panel of plots at the top in Fig. 2 shows the (pjj^ points and the fit where the fit 
function includes the lattice-spacing effects described in Sections 2 and 3. The plot at the bottom 
of Fig. 2 shows the extrapolations in the limit a = 0. The upper (D^) curve shows m/ m setting 
niq = nih = nis, while the lower (D+) curve shows mq,mi m setting nif, = nis. The physical quark 
mass inputs are from the MILC light meson analysis and m = {niu + md)/!. The points denoted by 
the red triangles correspond to physical /o, and Our preliminary results are: 

= 217 ± 10 MeV, fo, = 260 ± 10 MeV and /z),//d+ = 1.20 ±0.02 . (6.1) 

We have combined the statistical and the systematic uncertainties listed in Table 2 in quadra- 
ture. Our largest uncertainty is the combined uncertainty from statistical and residual discretization 
effects. The second largest uncertainty, chiral extrapolation, is an estimate of chiral expansion ef- 
fects not included in the fit function and effects from variation in the extrapolation procedure. The 
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mPiSq 



Figure 2: The preliminary D meson chiral extrapolation. The 3x4 matrix of plots (top) show the (j) data 
and corresponding fit including effects. Reading from left to right and top to bottom, plots correspond 
to {a,mi/mh) values of (0.15,0.2), (0.15,0.4), (0.15,0.6), (0.12,0.14), (0.12,0.2), (0.12,0.4), (0.12,0.6), 
(0.12,0.1), (0.09,0.2), (0.09,0.4) and (0.09,0.1). The larger plot (bottom) shows an overlay of the /o, 
and /£)+ extrapolations. The extrapolated curves are the fit (with error bands) taking ^ and fix- 
ing/extrapolating the light quarks to their physical masses. The extrapolations are not expected to go though 
any of the points which are computed at finite a. None of the data points having near seen the upper 
panel are visible in the Ds extrapolation. 6 
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Figure 3: Comparisons of the Fermilab/MILC values of /0+ and fo, to values from the HPQCD Collabora- 
tion [6] and recent experimental values [3] [5]. 



third largest error is the statistical error in the nonperturbative calculation of the current renormal- 
izations Zy and Zy^. The value of fo, is about eleven MeV (one sigma) higher than our earlier 
value. Using nominal kappa values rather than tuned values at the previous ri value accounts for 
about 1.3 MeV of the difference. Changing to the new ri while keeping kappa tuned results in a 4.2 
MeV increase. Incorporating heavy quark effects into the fit increases fo^ by about 2 MeV. Higher 
statistics on the most chiral of the a 0.09 fm lattice increases /o^ by about 1 MeV. These changes 
combine nonlinearly in the fit to yield the net increase. 

Figure 3 compares the Fermilab and MILC Collaboration values for the decay constants with 
the HPQCD Collaboration [6] values and with the experimental results. The experimental result 
for fj)+ is from CLEO [3] while the fo^ value is the Heavy Flavor Averaging Group average [5] 
of determinations by CLEO, BaBar and Belle. The Fermilab / MILC results remain in agreement 
with experiment. The total error on the experimental average for fj)^ is now smaller that our error 
providing a challenge for future lattice determinations. The apparent discrepancy between the 
HPQCD value of fo^ and the other two fo^ values is most striking. The HPQCD value is lower by 
about 1.8-2. la. The source of this difference may be clarified by further lattice simulations. 



7. Summary and future plans 

We have made several improvements in our analysis: a) discretization effects from both heavy 
and light quarks are modeled in our extrapolation function, b) we adopted a more precise ri value 
which derived from the MILC fji analysis rather than the ri value related to early T spectrum results 
c) we have improved the tuning of kappa charm. These improvements to the analysis will be more 
crucial in our next generation of decay constant study. We will increase statistics by a factor of four 
and extend the analysis to the finer lattice spacings a w 0.06 and 0.045 fm which will reduce our 
combined statistical plus discretization error as well as help reduce uncertainties attributed to chiral 
extrapolation procedures. In addition, a new high-statistics computation of the nonperturbative part 
of the current renormalization aims for an error below the 0.5% level. 
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